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In this talk we describe a procedure for isolating the infrared singularities present in gluonic 
scattering amplitudes at next-to-next-to-leading order. We use the antenna subtraction framework 
which has been successfully applied to the calculation of NNLO corrections to the 3-jet cross 
section and related event shape distributions in electron-positron annihilation. Here we consider 
processes with coloured particles in the initial state, and in particular two-jet production at hadron 
colliders such as the Large Hadron Collider (LHC). We construct a subtraction term that describes 
the single and double unresolved contributions from gluonic processes using antenna functions 
with initial state partons and show numerically that the subtraction term correctly approximates 
the matrix elements in the various single and double unresolved configurations. 
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1. Introduction 

In hadronic collisions, the most basic form of the strong interaction at short distances is the 
scattering of a coloured parton off another coloured parton. Experimentally, such scattering can be 
observed via the production of one or more jets of hadrons with with large transverse energy. In 
QCD, the scattering cross section has the perturbative expansion, 
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fi{x\)fj{x%)dx\dx<i (1.1) 



where the sum runs over the possible parton types i and j. The single-jet inclusive and di-jet cross 
sections have been studied at next-to-leading order (NLO) [|l|] and successfully compared with 
data from the TEVATRON. 

The theoretical prediction may be improved by including the next-to-next-to-leading order 
(NNLO) perturbative predictions. This has the effect of (a) reducing the renormalisation scale 
dependence and (b) improving the matching of the parton level theoretical jet algorithm with the 
hadron level experimental jet algorithm because the jet structure can be modeled by the presence of 
a third parton. The resulting theoretical uncertainty at NNLO is estimated to be at the few per-cent 
level [@]. 

In this talk, we will focus only on the NNLO contribution involving gluons and will drop the 
parton labels. At NNLO, there are three distinct contributions due to double real radiation radiation 



UO AWLO> 



, mixed real- virtual radiation da^ LO and double virtual radiation Aol,t LO , that are given 



V.2 



by 
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where the integration is over the appropriate Af -particle final state subject to the constraint that 
precisely m-jets are observed, 
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(1.3) 



As usual the individual contributions in the (ra + 2), (ra + 3) and (ra + 4)-parton channels are all 
separately infrared divergent although, after renormalisation and factorisation, their sum is finite. 
For processes with two partons in the initial state, the parton level cross sections are related to the 
interference of M-particle /-loop and 7-loop amplitudes [{^^ \^^')]m by 



d °NNLO ~ 
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In this talk, we specialise to the gluonic contributions to dijet production. Explicit expressions for 
the interference of the four-gluon tree-level and two-loop amplitudes is available in Refs. pi], while 
the self interference of the four-gluon one-loop amplitude is given in [|J]. The one-loop helicity 
amplitudes for the five gluon amplitude are given in [S]. This contribution contains explicit infrared 
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divergences coming from integrating over the loop momenta and implicit poles in the regions of 
the phase space where one of the final state partons becomes unresolved. This corresponds to the 
soft and collinear regions of the one-loop amplitude that were analyzed in The double real 
six-gluon matrix elements were derived in [^]. Here the singularities occur in the phase space 
regions corresponding to two gluons becoming simultaneously soft and/or collinear. The "double" 
unresolved behaviour is universal and was discussed in [g, g]. 



2. The antenna subtraction formalism 

There have been several approaches to build a general subtraction scheme for the double real 
contribution at NNLO [jl^, [11], [12|]. We will follow the antenna subtraction method which was 
derived in [jl3|] for NNLO processes involving only (massless) final state partons. This formalism 
is being extended to include processes with either one or two initial state partons [14, 15, |l6| ]. 

To render the contributions with different final states separately finite, the general structure of 
the subtraction terms at NNLO is 



d&NNLO - / {d&NNLO ~ d °NNLo) + / ddjW 



J NNLO 

P»i+2 -J09 m+ 2 
+ ( dd NNLO- dd NN 1 Lo)+ I d& NNW 

+ f *»$u» (2- 1 ) 

Jd<P,„ 

where da^ NL0 (do^^ ) is the subtraction term for the double radiation (real-virtual) contributions 
respectively. 

In this talk, we concentrate on the double unresolved subtraction term dd^ NLO relevant for the 
six-gluon contribution two-jet production in hadronic collisions. It is made up of several different 
contributions, that depend on how the unresolved partons are connected in colour space, 

d &NNLO = d &NNLO + d ®NNLO + d ®NNLO + d ®NNLO + d &NNLO 

(a) One unresolved parton but the experimental observable selects only m jets from the (m+ 1) 
partons. 

(b) Two colour-connected unresolved partons (colour-connected). 

(c) Two unresolved partons that are not colour connected but share a common radiator (almost 
colour-unconnected). 

(d) Two unresolved partons that are well separated from each other in the colour chain (colour- 
unconnected). 

(A) Correction for the over subtraction of large angle soft radiation. 

Contribution (a) is precisely the same subtraction term as used for the NLO (m+ 1 )-jet rate, 
and is the product of a three-parton antenna and reduced (ra+ l)-particle matrix elements. Con- 
tributions (b)-(d) are derived from the product of double unresolved factors and reduced (ra + 2)- 
parton matrix elements. Subtraction terms for these configurations can be constructed using either 



3 



NNLO real corrections to gluon scattering 



Joao Pires 



double soft limit for gg >! 

#PSpoints=IOOOO 
1487 outside the plot 

317 outside the plot 
59 outside the plot 




0.99997 0.99998 0.99999 1 
R 



1.00001 1.00002 1.00003 



(a) 

Figure 1: (a) Example configuration of a double soft event with sy : 
double soft phase space points. 



(b) 

' S\2 = s. (b) Distribution of R for 10000 



single four-parton antenna functions [ ]17| ] or products of three-parton antenna functions. For glu- 
onic processes, one encounters the four-gluon antenna F 4 ° for the first time. Finally, the large-angle 
soft subtraction terms (A) contains soft antenna functions which precisely cancel the remnant soft 



behaviour associated with the antenna phase space mappings for the final-final [18], initial-final 



and initial-initial configurations. Explicit formulae for the various contributions to dd^ NL0 are 



available in [19] 



Note that the subtraction terms are also needed in integrated form. When both radiators are 
in the final state, as needed for electron-positron annihilation, the integrated antennae are given in 



ref. [13]. For processes with one hard radiator in the initial state, the integrals are known [ 15 ] while 



the work is still in progress for processes with two hadronic initial radiators [16]. 



3. Numerical results 



To numerically test that the subtraction term correctly reproduces the same infrared behaviour 
as the matrix element, we generate a series of phase space points that approach a given double or 
single unresolved limit. For each generated point we compute the ratio, 



R 



dcr R 

uu awlo 

dd s 

NNLO 



(3.1) 



which should approach unity as we get closer to any singularity. 

As an example, we consider the double soft limit. A double soft configuration can be obtained 
by generating a four particle final state where one of the invariant masses sy of two final state 
particles takes nearly the full energy of the event s as illustrated in figure [j] (a). 

In figure [j](b) we generated 10000 random double soft phase space points and show the distri- 
bution of the ratio between the matrix element and the subtraction term. We show three different 
values of x = (s — Sij)/s [x = 10~ 4 (red), x = 10~ 5 (green), x = 10~ 6 (blue)] and we can see that for 
smaller values of x the distribution peaks more sharply around unity. For x = 10~ 6 we obtained an 
average of R = 0.9999994 and a standard deviation of a = 4.02 x 10~ 5 . Also in the plot we give 
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Figure 2: (a) Example configuration of a triple collinear event with s\jk — > 0. (b) Distribution of R for 10000 
triple collinear phase space points. 



the number of outlier points that lie outside the range of the histogram. As expected this number 
decreases as we approach the singular region. 

As a second example, we perform a similar analysis for the triple collinear limit with three 
hard particles sharing a collinear direction as shown in figure ^| (a). In this case the variable that 
controls the approach to the triple collinear region is x = s\jk/s. We show results for x = — 10 -7 
(red), x = — 10~ 8 (green), x = — 10~ 9 (blue). For 10000 phase space points with x = — 10~ 9 , we 
obtained an average value of R = 0.99954 and a standard deviation of a = 0.04. As before, the 
number of outliers systematically decreases as we approach the triple collinear limit. 

Similar behaviour is obtained for all of the remaining double unresolved limits. 



4. Conclusions 

In this contribution, we have discussed the application of the antenna subtraction formalism 
to construct the subtraction term relevant for the gluonic double real radiation contribution to dijet 
production. The subtraction term is constructed using four-parton and three-parton antennae. We 
showed that the subtraction term correctly describes the double unresolved limits of the gg — > gggg 
process. The construction of similar subtraction terms for processes involving quarks should in 
principle be straightforward. Together with the integrated forms of the antenna functions (see 
Refs. [|l5|] for the initial-final and Ref. [16] for the initial-initial configurations), these double real 



subtraction terms will provide a major step towards the NNLO evaluation of the dijet observables 
at hadron colliders. 
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